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Non-Oberbeck-Boussincsq (NOB) effects on the flow organization in two-dimensional 
Rayleigh-Benard turbulence are numerically analyzed. The working fluid is water. We 
focus on the temperature profiles, the center temperature, the Nusselt number, and on 
the analysis of the velocity field. Several velocity amplitudes (or Reynolds numbers) and 
several kinetic profiles are introduced and studied; these together describe the various fea- 
tures of the rather complex flow organization. The results are presented both as functions 
of the Raylcigh number Ra (with Ra up to 10 8 ) for fixed temperature difference A be- 
tween top and bottom plates and as functions of A ("non-Obcrbcck-Boussincsqness") for 
fixed Ra with A up to 60K. All results are consistent with the available experimental NOB 
data for the center temperature T c and the Nusselt number ratio Numob /Nuob (the la- 
bel OB meaning that the Obcrbeck-Boussincsq conditions are valid) . For the temperature 
profiles we find - due to plume emission from the boundary layers - increasing devia- 
tions from the extended Prandtl-Blasius boundary layer theory presented in (Ahlers et al. 
2006 J. Fluid Mech. 569, 409-445) with increasing Ra, while the center temperature it- 
self is surprisingly well predicted by that theory. For given non-Oberbeck-Boussinesqness 
A both the center temperature T c and the Nusselt number ratio Nunob/Nuob only 
weakly depend on Ra. 

Beyond Ra sa 10 6 the flow consists of a large diagonal center convection roll and two 
smaller rolls in the upper and lower corners, respectively ("corner flows"). In the NOB 
case the center convection roll is still characterized by only one velocity scale. In contrast, 
the top and bottom corner flows are then of different strengths, the bottom one being a 
factor 1.3 larger (for A = 40K) than the top one, due to the lower viscosity in the hotter 
bottom boundary layer. Under NOB conditions the enhanced lower corner flow as well 
as the enhanced center roll lead to an enhancement of the volume averaged energy based 
Reynolds number Re E = (^u 2 ) 1/2 L/v of about 4% to 5% for A = 60K. Moreover, 
we find Re-^Qg/ ReQ B « (/3(T c )//3(T m )) 1 / 2 , with f3 the thermal expansion coefficient 
and T m the arithmetic mean temperature between top and bottom plate temperatures. 
This corresponds to the ratio of the free fall velocities at the respective temperatures. 
By artificially switching off the temperature dependence of [3 in the numerics, the NOB 

f Present Address: Department of Mechanical Engineering, School of Engineering, The Uni- 
versity of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-8656, Japan. 

\ Present Address: Laboratoire de Physique de Ecole Normale Superieure de Lyon and CNRS, 
46 Allee d'ltalie, 69007 Lyon, France. 

% d.lohse@utwente.nl 



2 K. Sugiyama et al. 

modifications of Re E is less than 1% even at A = 60K, revealing the temperature de- 
pendence of the thermal expansion coefficient as the main origin of the NOB effects on 
the global Reynolds number in water. 



1. Introduction 

Rayleigh - Benar d (RB) conv e ction in the Oberbeck-Boussinesq (OB) approximation 
(|Oberbeckl (|l879h : iBoussinesql (|l903h ) is rather a mathematical concept than physical 
reality. On the one hand it is driven by a temperature difference A = Tj, — T t between the 
bottom and top plates, whose distance in height is L; on the other hand, the temperature 
dependences of the material properties such as the kinematic viscosity v, the thermal 
diffusivity k, the heat conductivity A, the isobaric specific heat capacity c p , and the 
isobaric thermal expansion coefficient f3 are all ignored in the OB approximation apart 
from that of the density p, which is assumed to vary linearly with T, 



p(T)np m (l-p m (T-T m )). 



(1.1) 



Here, T m = (T t + T{,)/2 is the arithmetic mean value of the warmer bottom plate tem- 
perature Tb and the colder top plate temperature T t , and p m and (3 m denote the density 
and the thermal expansion coefficient at the mean temperature T m , respectively. For- 
tunately, the OB approximation is rather good, if the material properties depend on 
temperature only weakly or if the temperature difference A between the bottom and 
top plates is kept small enough. However, if either the material properties are strongly 
temperature dependent (as e.g. the viscosity of glycerol in some temperature regime) 
or if the temperature difference between bottom and top plates is chosen to be large in 
order to achieve larger Rayleigh numbers, the deviations from the OB approximation 
are expected to become relevant. The c onsequences of these dev iations are called non- 
Oberbeck-Boussinesq (NOB) effects, see lAhlers et al. (|2006l . l2009h . To which extent they 
affect the Nusselt number and possibly could account fo r the differences between th e 
Oregon and t he Grenoble data is an ongo i ng controver s y Ichavanne et al. (fl997 , 2001 ) 



200.3); 



Roche et al] d200ll l2002h: INiemela et all (|2000| . l200lh : iNiemela k SreenivasanT 
Ashkenazi k Steinberd (jl999l )) 

The signatures of NOB effects studied in this article are (i) a deviation of the center (or 
bulk) temperature T c from the arithmetic mean temperature T m , (ii) a modified z-profilc 
of the area averaged temperature, which develops a bottom-top-asymmetry, (iii) different 
thermal boundary layer thicknesses ^ At at the bottom and top together with different 
temperature drops A& ^ A t across these BLs, (iv) a modification of the Nusselt number, 
best expressed by the ratio Nunob/Nuob, and (v) a bottom-top symmetry broken flow 
structure, in particular possibly different amplitudes Ub ^ Ut of the plate-parallel winds 
near the bottom and top plates meaning different Reynolds numbers Reb,t = Ub.t/vL^ 1 . 

The deviation T c — T m of the center temperature from the arithmitic mean temperature 
and the corresponding differences between the temperature drops A t and A& over the 
thermal BLs presumably i s the e xperime ntally most explored NOB effect, nam ely for 
water by IWu fc Libchaberl (Il99ll) and by lAhlers et all (120061 . 120071 . 12008L l2009h . While 
for water even for A = 50 K th e deviation T r — is at most 2 K , for g lycerol this 
deviation can be as large as 8 K ( Zhang et al. (1997), Sugivama et all ( 2007f )). 



Ahlers et all (|2006l ) the Prandtl 



To theoretically account for these deviations, in 
Blasius BL theory was extended to the NOB case, giving surprisingly good agreement 
with the experimental center temperatures T c for water. Also for ethane gas, which is 
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compressible, the BL theory - extended to compressible fluid flow - can describe t h e mea - 
sured center tempera ture data rather sat isfactorily, as presented in lAhlers et all (|2007l ). 



But as was shown in AM ers et all (|2008l ) , the T-dependence of the buoyancy caused by 



thermal expansion (3 = j3(T) is the dominant cause of the observed NOB effects, in par- 
ticular the characteristic differences between the more gas-like and the more liquid-like 
ethane on the two sides of the critical isochore. 

The success of the Prandtl-Blasius BL theory in the context of NOB convection is 
remarkable for at least three reasons: First, the boundary layer theory deals with semi- 
infinite plates, while experiments are done in finite aspect ratio containers, mainly for 
r = 1. - Second, and more importantly, the Prandtl-Blasius BL theory completely 
ignores the plume separations and the corresponding time dependence of the boundary 
layer flow. Although the shear Reynolds numbers in the BLs are not yet very large in RB 
flow (for Pr = 1 the transitional shear Reynolds number Ref « 42 0, which indicates the 



range of turbulence transition, is only reached near Ra ~ 10 , see iGrossmann fc Lohse 



( 2002f )). the plumes ( and thus the time depe n dence of the BL flow) play a significant role 



in the heat transfer ( Ciliberto et al. ( 19961 ): Ciliberto fc Laroche ( 19991 )) and perhaps 



also for the bulk temperature T c . Plumes are not included in the classical BL theory 
because that does not take notice of the buoyant forcing in the Prandtl approximation of 
the hydrodynamic equations of motion. Neither the thermal expansion coefficient f3 itself 
is addressed nor is its temperature dependence taken into account, which in reality is 
considerable. For water at T m = 40°C and A = 40K there is nearly a factor of 2 between 
the respective values j3 t and /3j, at the top and bottom plates. Technically speaking, the 
BL theory misses (3 since buoyancy shows up in the equation of the vertical velocity field 
u z ; and this equation enters into the Prandtl approximation only to derive the height 
independence of the pressure, which in RB anyhow does not play a role. The numerical 
simulations presented in this paper will quantitatively show that buoyancy and plumes 
indeed affect the temperature BL profiles. 



The re are two further assumptions of the extended BL theory developed in lAhlers et al 



( 20061 ) which need to be tested: First, the extended BL theory assumes that the large scale 



wind velocity is the same close to the top and the bottom plate, i.e., Ut = Vb or Ret = 
Reb = ReNOB- We will show that there are several relevant velocity amplitudes, most of 
which break the top-bottom symmetry. Only the main central roll has the same amplitude 
near the bottom as well as the top boundary layers. Second, within the BL theory the 
ratios Nunob/Nuob and Re^oB I ReoB cannot be calculated at all, in contrast to T c . 
These quantities can only be obtained by additional input from experiment, namely by 
employing the experimental information on the ratio F\ := 2X S Q B /(Xf + Xf). Here A*' 
and A^ are the top and bottom thermal BL thicknesses, defined via the tem perature 



slopes at the plates in Eq. (|4.ip and sketched in figure [TJ In lAhlers et all (|2006l ) we have 
calculated F\ from the measured Nusselt number and the calculated center temperature. 
Its value T c determines the ratio Fa := («tA( + K;,A&)/(k to A)). The exact relation holds 
(even for compressible flowjj] 

Nu NO b 2A^ B K t A t + K b Ab 
Nuob Af + Xf K m A 

Nuno b is the actual heat flux with all material parameters taken at their respective real 
temperature values. The label OB means that all fluid properties are taken as temperature 
independent constants, evaluated at the arithmetic mean temperat ure T m . 

Remarkably, for the analyzed case of water at T m = 40°C (see Ahlers et al. ( 20061 )) 



t The notation used in eq. (|1 .2|) is explained in the caption of figure [T] 



4 



K. Sugiyama et al. 



the ex perimental data are consistent with Fx = 1. But as was shown in lSugivama et 

by numerical simulations it is Fa ^ 1 for RB convection in glycerol under NOB 
conditions, at least up to Ra = 10 8 . However, in glycerol due to the large Prandtl number 
Pr = 2500 a large scale convection roll did not yet develop and it could be that Fx = 1 
is connected with the existence of such a roll. The numerical simulations presented in 
this paper will unambiguously show that F\ = 1 does not hold in general. This property 
thus is co-incidental for water due to the specific temperature dependences of its material 
properties around 40°C. 

The question of modifications of the Reynolds number(s) through NOB effects is inti- 
mately related to the flow organization. In recent years there was consider able progress in 
our insi ght into the flow struc t ure, thanks to nume ri cal simulations (see e.g . Verzicco &: Camussi 



(l2003l); ISchmalzl et all (|2004|) ; I Amati et all (l2005l ); IStringano fc VerzicccT(l2006l)L to PIV 



measurements (|Xia et all (|2003l ); lQiu et all (l2004D;IXi et all (l2004l ); ISun et all (|2005aibh ). 



and to velocity correlation measurements ( Brown et all ( 2007t )). These papers revealed 
that there are various feasable possibilities to define flow amplitudes and that these 
differently defined amplitudes and the corresponding Reynolds numbers have different 
scaling behavior with Ra. Our numerical simulations have fully confirmed and detailed 
this view. We will show that NOB conditions influence the flow structure near top and 
bottom differently and modify the various Reynolds numbers correspondingly. NOB 
conditions have the largest impact on the convective flow in one top and one bottom 
corner of the cell, where macroscopically visible secondary rolls develop. We will also 
show that the NOB modification of the global, volume and time averaged, energy based 

1 /2 

Reynolds number Re E = {\ u2 ) Vt l( v L~ v ) is consistent with attributing it mainly to 
the change of the thermal expansion coefficient (3 in the bulk. More specifically, we find 
R e NOB/ R e OB ~ (/?(Tc)//3(Tm)) 1 ^ 2 , a finding clearly not describable within the extended 
Prandtl-Blasius BL theory. 

In this paper we focus on water (Pr — 4.4). Nevertheless, the parameter space is 
considerable. Next to Ra the crucial parameter is the NOBncss A. For comparison we 
perform numerical simulations for fluids with non-physical temperature dependences of 
their material properties in order to clarify the origin of certain observations. 

As the numerical effort is so large for three-dimensional simulations we restrict our- 
selves to two-dimensional simulations. One may worry on whether two-dimensional simu- 
lations arc sufficient to reflect the dynamics of three-dimensional R B convection. F o r heat 
flow under OB conditions thi s point has been analyzed in detail by Schmalzl et all ( 20041) 
and e arlier by iDeLuca et al. (ll990h : IWerne et ali (Il99lh : IWernel (09931 ). ISchmalzl et al 



( 20041 ) 's conclusion is that for Pr > 1 various properties observed in numerical 3D con- 
vection (and thus also in experiment) are well reflected in the 2D simulations. This in 
particular holds for the BL profiles and for the Nusselt numbers. Since one focus of this 
paper is on the difference between OB and NOB convection, the restriction to 2D simula- 
tions might be even less severe, as NOB deviations occur in both cases and the differences 
between 2D and 3D simulations might cancel out in quantities such as (T c — T m )/T m , 
Nunob/Nuob, or Re nob I ReoB- We also note that for a comparison with the Prandtl- 
Blasius BL theory 2D simulations are in fact more appropriate than 3D simulations, as 
the BL theory is two-dimensional per construction. 

The paper is organized as follows: In section [5] we will explain, justify, and verify 
the numerical method. Section [3] is devoted to our results on the mean temperature 
profiles and the related shifts of the center temperatures. Section [4] addresses the NOB 
effects on the Nusselt number. The main section is section [5] where we first analyze 
the flow structure for the OB case and then its modifications through NOB effects. 
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z = L 



Figure 1. Notations and sketch of the time and surface averaged temperature profiles versus 
height z in the OB and NOB cases, respectively. The height of the cell is L (not shown). The 
temperatures at the top (z = L) and the bottom (z — 0) plates are T t and T),, their arithmetic 
mean is T m = (Tt + Tt)/2. The thermal BL thicknesses based on the temperature slopes at 
the plates are called \f b . The respective temperature drops are At and A&. The time averaged 
temperature in the bulk (or center) is T c . For water as the working fluid this center temperature 
T c is larger than T m . While A^ = \f in the OB case, under NOB conditions the bottom BL is 
thinner than the top one, \ f < Af . The fluid properties such as u, k, and (3 carry the same index 
as the temperature at which they are considered, e.g. v t = v{Tt) for the kinematic viscosity at 
the top plate, and so on. 



Several feasable measures for the wind amplitudes of the complex flow structure will be 
introduced. Section [5] contains the conclusions. 



2. Definitions, governing equations, and numerical method 

The equations governing non-Obcrbeck-Boussincsq convection in incompressible liq- 
uids are the incomprcssibility condition 



dim = , 



(2.1) 
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Table 1. Expansion co efficients of materia l properties of water around the temperature 
= 40° C adapted from |AM ers et a/1 (|2006j) . The kinematic viscosity, the thermal diffusiv- 
ity, and buoyancy are written in a polynomial form as u(T) = t}{T) / p m ~Yl =o A n {T — T m ) n . 
k(T) = A(T)/(p mCp , m ) = J2 n=0 B n (T - T m ) n , and g(l - p(T) / p m )=Z n =o C «( T " T ™)"> re " 
spectively. Using the leading coefficient Ci(= /3 m ) for the buoyancy force, we can write the 
Rayleigh number defined in (|2.4|) as Ra = CiL 3 A/(v m K m ), where v m = rj{T m )/pm and 
K m — A(T m ) / (pmCpm) , which coincides with the usual OB definition. 



the Navier-Stokes equation 



p m (d t Ui + UjdjUi) = -dip + dj(rj(djUi + diUj)) + p m g (1 — — ) 5 i3 , 

Pn 



and the heat-transfer equation 



j d j T) = d j {kd j T). 



(2.2) 



(2.3) 



The dynamic viscosity rj(T) and the heat conductivity A(T) are both temperature and 
thus space dependent. The isobaric specific heat capacity c p and the density p are as- 
sumed as constants and their values c p ^ mi p m fixed at the temperature T m , except in the 
buoyancy term, where the full nonlinear temperature dependence of p(T) is implemented. 
For water, which we here consider as the working fluid, density as well as specific heat 
are indeed constant with temperature to a very good approximation. The experimentally 
known temperature dependences of 77, A, and p (in the buoyancy term) t ogether with 
the va lues of the parameters p m , c p ^ m for water are given in the appendix of lAhlers et all 
and, for better reference, are reported in Tab. [T] in the form implemented in the 
present DNS. 

We deal with a wall-bounded system with an aspect ratio fixed at T = 1. The veloc- 
ity boundary conditions accompanying the dynamical equations are Ui = at the top 
and bottom plates z = L and z = as well as on the side walls x = 0,x = L. The 
temperature boundary conditions are Xj, — T t = A for the temperature drop across the 
whole cell of height L. At the side-walls (x = 0, x = L) heat-insulating conditions are 
employed, d x T\ x= Q^ = 0. The cell is considered to be 2-dimcnsional, i.e., there is no 
y-dependence. The Rayleigh number is defined with the material parameters taken at 
the mean temperature T m , 

Ra=^ L3A . (2.4) 

We vary the Rayleigh number in DNS by varying the height L of the box, while the non- 
Obcrbcck-Boussincsqncss is changed by varying the temperature drop A. Note that in eq. 
(|2.2p the full temperature dependence of the density in the buoyancy term is taken into 
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account, rather than employing the linear approximation eq. (jl.ip only. Still the Rayleigh 
number is defined with the coefficient of the linear expansion of the density with respect 
to temperature, taken at the mean temperature, /3 m — p^dp/dT\T m , cf. Tab. [1] The 
Prandtl number Pr = v m / ' n m is also defined in terms of the material parameters at the 
arithmetic mean temperature. 

Equations (|2.ip - (|2.3[) are solved on a two-dimensional domain with gravity pointing 
in negative ^-direction. To discretize the N avier-Stokes and heat t r ansfer equations, 
we em ploy a finite difference scheme (see e.g. Pevret fe Taylor (1983|); Ferziger fc Peric 



1 19961) 1 The space derivative s are approximat e d by the fourth-order central difference 



scheme on a staggered grid ( Harlow fc Welch d 19651)). In particul ar for the advection 



terms we employ the scheme proposed by iKaiishima et al. ( 2001 ). which satisfies the 



relations dj(UjUi) — UidjUj + UjdjUi and dj(ujT) — TdjUj + UjdjT in a discrctized 
form and ensures that the second moments of the velocity and temperature are highly 
conserved. To integrate the equations in time, we use the second-order scheme, i.e., the 
Adams-Bashforth method for the advection te rms and the C r ank-N icolson one for the 
viscous, diffusive, and buoyant terms (see e.g. ICanuto et all ( 19881 )) . To complete the 



time marching in the momentum equation and simultaneously satisfying the solenoidal 
condition (12.11) of the velo city vector, we employ a simplified-marker-and-cell procedure 



( Amsden fc Harlowl ( 1970l )) by solving a Poisson equation for the pressure. The two- 



dimensional discretized pressure equation, which is written in the fourth-order finite 
difference form, is reduced into a one-dimensional problem by taking the Fast Fourier 
Transform (FFT) in the x-direction. The boundary condition at the side walls (x = and 
x = L) is satisfied, if the relation d x 4>\ x =o = dx4>\x=L = holds for all the quantities (f> 
in the pressure equation. To impose this condition, we take a periodicity 2L for the FFT 
(i.e., 4>{x) = <p(2nL + x) with n an arbitrary integer) and introduce a fictitious domain 
L < x < 2L, in which the quantities are given by <p(x) — cf>(2L — x). We directly solve 
the reduced-order equation written in a heptadiagonal matrix form and then determine 
the pressure field by taking the inverse FFT. 

We have validated the numerical code by checking the instantaneous kinetic energy 
and entropy budget relations for u 2 /2 and T 2 /2 both in the OB and NOB cases and by 
evaluating the correctness of the onset of convection in the OB case. The critical Rayleigh 
nu mber we compute ( Ra r = 2585.27) is in agreement with the one computed analytically 



bv lLuijkx fc Plattenl (|l98ll ) (Ra c = 2585.03) to a precision of less than 0.01%. We note 



that Ra c is much larger than the more known Ra c = 1708 for an infinite aspect ratio 
system, due to the presence of lateral walls. 

The area averaged heat currents are calculated as functions of time r both at the top 
(t) and the bottom (6) plates separately, 

-Ktd z (T) At (z = Wr) ~K b d z (T) Ab (z = 0,r) 

^ (r) = ^aTZ ' Nub{T) ^aTZ ' (2 - 5) 

with (. . .}a denoting the averaging over the horizontal surface (actually over x only, since 
the system is 2D) of the top and bottom plates. We find good agreement of the time 
averages of Nu t (r), Nub(r), see figure [2l 

To quantify the statistical convergence, we make an uncerta i nty an alysis estimating 
the time autocorrelations (see for instance Tennekes fc Lumlev ( 19721) Sect. 6.4). For a 



time-dependent function /(r), the error is evaluated as Sf = / rms yj 2ti / T ; here / rms 
is the root mean square of /, 77 the integral time obtained from the autocorrelation 
coefficient of /, and T the total simulation time under statistically steady conditions, 
i.e., after a transient time tq. Considering the error propagation, we evaluate the errors 
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Figure 2. Upper two panels: Temporal evolution of Nut(r) and Nui,(t) as defined in equation 
(|2.5[) . before time averaging, for Ra = 10 7 , A = 40K, and L = 1.89 cm. The total integration 
time T, reported here in seconds, is of the order 10 4 integral times of the autocorrelation 
coefficient, ri ~ 0.606 s. Lowest panel: Ratio of the temporally averaged Nusselt numbers 
f£ Nut(r) dr/ J* Nub(r) dr. Here the time averages at top and bottom agree up to 0.009%. 
For other Ra and A the statistical convergence is similar with maximum relative error being 
0.091%. 



of F\ ■ Fa, of F\, and of Fa at 95% confidence level as, respectively, 



S(F X F A ) - Nunob ( 6 ( Nu nob) \ | f S(Nu OB ) \- 
Nuqb V V Nunob J \ Nuqb ) 



S(F- 



)=p ( S(Nu NOB ) y ( S(Nu OB ) y (S(F A )\ 
' A V V Nunob J \ Nu OB J \ F A J 



S(F A ) 



S(T C ), 



(2.6) 



which will be indicated by the error bars in the plots. 
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Figure 3. Temperature profiles 9(£) = ({T} A tf) - T t )/A for Ra = 10 4 and A = 40K for water 
at T m — 40°C. The symbols indicate the results of the numerical s imulation, th e lines stem 
from the extended Prandtl-Blasius boundary layer theory proposed in I Aiders" et aZ.I(|2006r ). The 
numerical simulations give O(z). £ = z/Ibl is the similarity variable. In classical BL theory its 
scale is given by Zbl = \/xu m /Ut,b with x being the distance from the plate's edge; in this figure 
in order to translate to 0(£) we have chosen the factor Ibl such that the curves have the same 
slope at f = . 



3. Mean temperature profiles and center temperature 

We will first focus o n the water case with T m = 40°C, corresponding to the experiments 
of Ahlers et al. (|2006h . In figures [3] and g] the mean temperature profiles, averaged over 
the full width < x < L of the cell, are shown for Ra = 10 4 and Ra = 10 8 , respectively 



and c ompared with the extended Prandtl-Blasius BL theory developed in lAhlers et 



( 20061 ) (the specific procedure adopted for the comparison is detailed in the caption of 
Fig.©. 

For Ra = 10 the agreement between numerical data and the time independent ex- 
tended BL theory is excellent, both for the profiles and for the center temperature. This 
good agreement is remarkable as originally the Prandtl-Blasius BL theory has been de- 
rived for semi-infinite or at least long flat plates. 

For Ra = 10 s there are differences between the numerical and the BL theory profiles, 
namely, the numerical profiles are somewhat smoother than those from the BL theory, 
upper panel (a) of figure 3J At this relatively large Ra such differences are not unex- 
pected, because of the enhanced plume activity at larger Ra, which is not included in the 
extended BL theory. The plume activity is not homogeneous in the horizontal direction 
and the area averages in evaluating 0(£) in the upper figure is taken over the whole 
container width < x < L, including also the neighborhoods of the side walls, where 
the plume convection is preferentially strong. The influence of the plume convection can 
be confirmed by comparing with the numerical, only time averaged z-profiles along the 
middle line x — L/2, see lower panel of figured! These profiles now show good agreement 
with the extended BL theory , much better than those in the upper panel. We attribute 
this to the expectation that the plume detachment near x — L/2 is less than that near the 
side walls, and thus the BL approximation should be more reasonable along the middle 
line. It could be objected that the observed deviations might be due to the influence of 
the side- walls, which are included in the horizontal surface averages but excluded for the 
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(a) 



1.0 




Figure 4. Temperature profiles at Ra = 10 with A = 40K and T m = 40°C in water. Upper 
panel (a): horizontally area averaged temperature profiles as in figure [3] Though the center 
temperature T c is well described by the extended Prandtl-Blasius BL theory, for this relatively 
large Ra the numerical temperature ^-profiles show significant deviations from those of the 
extended BL theory. We attribute these deviations to the plume detachments, which are not 
included in the extended Prandtl-Blasius BL theory. To demonstrate this we show in the lower 
panel (b): time averaged temperature z-profiles at a fixed x- value, namely along the middle line 
x = L/2, i.e., 0(£)U=l/2 = {T\ x= l/2 — T t )/ A. Here the plume activity is expected to be weaker 
than in the regions near the side walls, which contribute to the area averaged profile in panel 
(a). Indeed DNS and BL theory agree satisfactorily along the center line. 



center-line, with only time averaging. To clarify this point we have performed a DNS with 
lateral periodic boundary conditions and T = 2. The numerically obtained profiles (not 
shown here) are very close to the surface and time averaged profiles of Fig. 0Ja) rather 
than to the time averaged center line profiles in Fig. 0Jb) . This supports the conclusion 
that it is the plume flow and temporal dynamics of the BLs, which is the main reason for 
the observed discrepancy between the Prandtl-Blasius BL theory and the area averaged 
DNS profiles. 

Nevertheless, in spite of the large deviations in the z-(or £-)depcndence of the area 
averaged temperature profiles as shown in figurelUa), the center temperature T c obtained 
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Figure 5. Relative deviation (T c — T m )/A of the horizontally area (and time) averaged center 
temperature T c from the arithmetic mean temperature T m in terms of A versus the temperature 
difference A for water at fixed T m = 40°C for various values of Ra. The nearly linear increase 
means that T c — T m = constA 2 + h.o.t. 
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Figure 6. (T c - T m )/A vs Ra for water at fixed T m = 40°C and A = 40K. Inset: deviation of 
T c from the center temperature in the laminar case Tc" 1 ' 1 at the onset Ra c = 2590 of thermal 
convection at T = 1. Note the strong variation of T c in a range slightly above Ra = 10 , which 
we attribute to transitions between different coherent flow structures. 



from the extended BL theory ( Ahlers et al. ( 200ot )) still very nicely agrees with that 
calculated with DNS. This is confirmed in figure [SI showing T c as a function of the 
NOBncss A for various Ra, ranging from 2 ■ 10 3 to 10 8 . In figure [5] we display the T c - 
shift for fixed A = 40K and T m = 40°C as a function of Ra. Interestingly enough, 
beyond some 10 5 for Ra, the center temperature T c is rather independent of Ra. Only in 
the immediate range beyond the onset of convection. T c — T m is pronouncedly smaller, 
reflecting the smooth transition to the small value of T c — T m in the non-convecting state. 
The Ra dependence of T c is not monotonous, what we attribute to transitions between 
the various coherent RB flow patterns in the considered i?a-range. 
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Figure 7. Nusselt number ratio Nunob /Nuob vs. Rayleigh number Ra for water at fixed 
values for T m = 40°C and A = 40K. It is also indicated where we are in the steady, periodic, or 
chaotic regime. 



The results shown in Fig. [6] arc consistent with previous findings (ILohse fc Grossmann 
( 19931 )) that in the i?a-range from onset of convection up to about Ra « 5 • 10 7 the flow 
only successively looses its spatial coherence. In this i?a-range the relative coherence 
le ngth (-coherence! L decre ases from values far above 1 to values of order 1/6, see figure 1 
of Sugiyama et al\ ( 2007t ) in wh ich we have calculated (coherence i 1 L as a funct ion of Ra, 
based on the unifying theory of iGrossmann fc Lohsel (|2000l l200li 120021 . l2004h . Only for 
Ra beyond this transition range up to some 10 7 , in which spatially coherent structures 
are gradually lost, the heat convection is fully turbulent and the pdf of the fluctuations 
becomes exponential instead of being Gaussian. 



4. Mean heat flux 

Similarly to T c , the Nusselt number ratio Nunob/Nuob displays an only weak de- 
pendence on Ra, see figure [7] This ratio has been written in (|1.2p in terms of two factors, 
Fa and F\. The latter one measures the changes of the thermal BL widths caused by 
NOB conditions. These widths of the temperature profile are sketched in figure [TJ quan- 
titatively we define them in terms of the temperature slopes at the plates. 

\ sl = At \ sl = Ab (A i ^ 

' ~\d z (T) AuT y Ab ~ \d z (T) Ab ,r\- [ ■ ' 

In figures [S] and [5] we reveal the origin of the Nusselt number modification in the NOB 
case. In general it is a combination of the A-dcpcndcncc of both factors F\ and F& 
in the product (|1.2p . As shown in figure [5] for fixed non-Obcrbcck-Boussincsqness A, 
Fa displays a weak dependence on Ra for Ra > 10 4 . This can be understood from the 
weak i?a-dependence of the center temperature T c on the NOB changes of the material 
parameters (see figures [5] and [6]) , and because Fa = ((k* — Kb)T c — KtTt + KbTb) / (K m A) 
depends on T c only (for given Tt and Tb). 

On the other hand the factor Fx, describing the variation of the thermal BL thicknesses, 
shows a rather weak but obvious dependence on the RB flow regimes. For the fully chaotic 
regime (Ra > 10 6 ), the deviation of Fx from Fa = 1 is much smaller than that of Fa, 



Flow organization in Rayleigh- Benard convection 



13 



1.010 



1.005 



Periodic 



^ 1.000 • 
0.995 



\ Steady 

X x xxxxX xX 



-X— ' Fx 



0.990 



4*> 



BL theory F A 
I 

^AA- A- --^ - -A - A- - 
_■ i_i ■ ■ i_i i ■ i_i 

10 J 10 4 10 5 10 6 10 7 10 8 
Ra 



Chaotic 



11 L H 



Figure 8. F\ and Fa (defined in (|1.2|l and in the text before) versus Rayleigh number Ra for 
water at fixed values for T m = 40°C and N OBness A = 40K. T he dashed line corresponds to Fa 
resulting from the BL theory developed in lAhlers et al\ (|2006f ). Reasonable agreement between 
BL theory and DNS is observed. Fx in the presently considered water case is compatible with 
Fx = 1 for larger Ra, as in experiment. The significant Fa-dependence for smaller and medium 
Ra seems to reflect the changes of the coherent flow structures still present at these Ra. In 
particular one recovers the window of quite different behavior slightly above Ra = 10 . 
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Figure 9. Nusselt number ratio Nunob / Nuob = F\ - Fa together with the factors Fa and Fa 
versus NOBness A for fixed Rayleigh numbers (a) Ra = 10 4 , (b) Ra = 10 7 , and (c) Ra = 10 8 . 
The working liquid is water at T m = 40° C. Th e dashed lines correspond to Fa resulting from 
the BL theory developed in lAhlers et al.\ (|2006f ) . 



which indicates that Nunob/Nuob is dominated here by Fa and thus by the behavior 
of the center temperature. As shown in figure [51 which displays the dependences on 
the NOBness A for the largest analyzed Ra = 10 7 and 10 8 , the A-dcpcndence of F\ 
happens to be very small, and for those Rayleigh numbers F\ « 1 happens to be a 
good approximation. This might be due to an incidental combination of the temperature 
dependences of the material parameters tj(T), A(T) and p{T) around the chosen mean 
temperature T rn = 40°C in the case of water. The experimental finding, reported in 
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Figure 10. Same as figure |9j but now with a temperature independent thermal expansion 
coefficient /3m, i- e., with only linear temperature dependence of p(T) as in eq. (|l.l|l . The 
A-dependences of the factors Fa and Fx and thus Fa ■ Fx are different from those in the 
case of fully temperature dependent P{T). 



Ahlers et al. ( 20061 ). that F\ ~ 1 in a similar Fa-range, for the same T m = 40°C, and 
for A up to 40K therefore can be considered as incidental. It is not a general property 
of NOB Raylcigh-Bcnard convection. We have also checked the influence of buoyancy on 
the BL widths as described by F\. If we disregard the T-dependence of /3(T), i. e. the 
non-linear temperature dependence of p(T) in our numerical simulations and thus have 
a constant (3 — /3 m , F\ shows a larger deviation from 1 at Ra = 10 and happens to be 
closer to 1 at Ra = 10 4 , sec figurcQIJl just opposite to the case with full T-dcpendence of 
/3(F). Still the nonlinear temperature dependence of p, i.e., the temperature dependence 
of P, has a relatively weak effect on F\ and Fa and therefore on the Nusselt-number 
modification. 

The second conclusion we can draw from figures [8 ) [9) and [TOl is tha t in all cases Fa 
is correctly described by the BL theory given in Ahle rs et aZT(|2006[ ). As F A can be 
calculated from F c only, this is of course to be expected, since T c is well described by the 
extended Prandtl-Blasius BL theory (see figure [5|) . 



5. The flow structure and various wind amplitudes 

We focus now on the structure of the flow field or "wind" in thermal convection. 
Although in our 2D simulati ons we miss interesting but typically 3-dimensional flow 
modes (cf. lAhlers et al. ( 20091 ) for a summary), even in two dimensions the velocity field 
is rather complex, as the snapshots in figure [TT] display. Nevertheless such convection 
fields have mostly been described by only one single a mplitude U. This wind amplitude 



U is a cruci al parameter both in the general OB theory (jGrossmann fc Lohse (I2000LI2001 



2002l|2004l) ) as well as in the extended BL theory dealing with NOB effects ([Ahlers et al 



(|200fii )T The GL theory hypothesizes that one needs only one single mean wind amplitude 
to describe the heat transport (also for large Ra) and that this amplitude is essentially 
uniform throughout the cell. 

The extended BL theory, developed in Ahlers et al. ( 20061 ) for NOB situations, as- 
sumes that such a uniform wind is still present even under NOB conditions and that in 
particular the top and bottom BLs see the very same wind amplitude Ut = Ub = Unob, 
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Figure 1 1 . (color) Snapshots of the velocity (arrows) and temperature (color) fields for Ra — 1 8 
at T m — 40°C, working fluid water. The left panel corresponds to the OB case (all material 
properties are kept temperature independent, taken at T m ), the right one to the NOB case, 
both with the same A = 40K. The temperature color scheme is in °C, same in both panels. 



in spite of the NOBness. The amplitude Unob is allowed to be different from Uob ■ but 
for the Nusselt number calculations (not for T c , as detailed above) its value has to be 
taken as a parameter of the theory. It is this character of Unob as a boundary condition 
for the Prandtl-Blasius BL equations which leaves the BL theory incomplete for calcu- 
lating the heat transport across the RB cell. Thus neither the Reynolds nor the Nusselt 
number deviations under NOB conditions can be predicted, unless further input (data 
or assumptions) is introduced (as for instance F\ = 1). 

We organize our analysis of the 2-dimcnsional flow structure as follows. Starting with 
the visualization of the dynamical flow fields, we next introduce a time averaged con- 
vective Eulerian field. Then we discuss several sensible possibilities to adequately define 
relevant wind amplitudes quantitatively. Finally, we present our results about the depen- 
dences of the [/-amplitudes or Re = U / (vL^ 1 ) numbers on the Raylcigh number Ra and 
the NOBness A. 



5.1. Dynamical features 

For RB convection in water in an aspect ratio Y = 1 container it is known from experiment 
and also noticeable in our 2D DNS flow as shown in the snapshots in FigfTTl that a large 
scale circulation (LSC) with an extension comparable to the box size L is present both 
in the OB and NOB cases at Ra = 10 s . In addition and on top of this LSC there are of 
course fluctuations of the u, T-fields. 

The time development of the numerical ta-fiel d shows reversals in the c ir culation sign 



in agr e ement with experimental ob s ervations dSreenivasan et al\ ()2002t ) : iBrown et al 
(|2005l);lBrown fc Ahlersl (l200d 12007^1 Ahlers et all (|2009l) ) and earlier numerical simula- 
tions ([Hansen et al\ ( 1992h ). In 3-dimensional experiment these reversals can occur either 
by rotation of the convection roll's plane or by cessation and restart; in the 2-dimensional 
numerics of course only the latter type of reversal occurs. Several models have been devel- 
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oped for the s e reversals, see e.. 



Sreenivasan et all (|2002l ); lFontenele Arauio et all (|2005l ): 



Benzil (|2005h : iBrown fc AhlersI ^OOTt ). Wc shall report on details of our results about the 



statistical properties of the reversals in our 2D simulations elsewhere. In the context of 
the present paper these reversals only complicate the statistical analysis of the flow field, 
as long-time averages of the velocity field become zero and wash out the flow structures. 

How then to obtain the main features of the dynamical, complex time-dependent u- 
field? To achieve them we consider conditionally time averaged Eulerian fields as well as 
several profiles, which partly take the fluctuations into account too. In addition a global, 
energy based wind amplitude is introduced. 

5.2. Conditionally time averaged velocity fields 

To overcome the problem that long time averages due to the statistical flow reversals 
give zero velocity everywhere, we perform conditioned time averages, which take the time 
dependent rotational direction of the wind into account. This instantaneous rotational 
direction is identified by the sign of the vorticity at the center of the box. Whenever the 
wind reverts its direction, before performing the standard time-averaging the velocity 
field is mirrored along the vertical center-line. Respecting this, from the full velocity field 
u(x, t) we can compute another, time averaged, complete 2D Eulerian-type velocity field 
u(x), in which the respective local direction of the velocity is coupled to the sign of the 
central vortex. Component-wise we define this conditionally time averaged flow field as 

j-to+T 

u^(x,z) = -T _1 / dt u x (x,z,t) sign uj c (t), 
u z {x,z) = T 1 / dtu z (x,z,t). 

J t 

Here T denotes the averaging time, sign u) c (t) is the sign of the vorticity u)(x, z) = 
d z u x — d x ii z at the center of the box (x, z) = (4, -|), and x is 

x = L/2 + (L/2 - x) ■ sign uo c (t) . 

In a similar manner we also define the conditionally averaged velocity squares as 

i-ta+T 

u x (x,z)=T~ 1 dtul(x,z,t), 

_ \ + r (5-2) 
v%(x, z) = T 1 / dtul(x,z,t). 
J t 

Rather than conditioning on the sign of the vorticity at the center, on e could also con 



dition on the sign of the total angular momentum, as has been done by Ivan Heiist et 



(l20f)6h (which minimizes the contributions of some high-frequency oscillations), but for 



the purpose of this paper the difference between these two types of conditional averaging 
have turned out not to be relevant. 

Conditionally time averaged fields u(x), being time independent objects, allow to visu- 
alize the persistent spatial structures in the flow field. Some conditionally time averaged 
fields obtained for different Ra numbers are shown in FigfT2l 

Note again the large circulation roll in the center range but also the secondary counter- 
rotating rolls in the corners. We remark that suc h secondary circu lation rolls in ncarly-2D 
convection have been experimentally detected bv lXia et all (|2003l ). In our simulations sec- 



ondary rolls appear for Ra > 10 5 . We interpret the secondary rolls as caused by boundary 
layer separations, which are known to occur when a flow is heading a perpendicular wall. 
They might be considered as kind of "wakes" behind the separation. Remarkably, the 
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Figure 12. Lines of constant values for the conditionally time averaged velocity field u(x,z) 
at different Rayleigh numbers Ra = 10 4 , 10 6 , 10 s in OB convection. Counterclockwise veloc- 
ity direction is drawn with solid lines, while the clockwise ones are indicated by dotted lines. 
The OB flow structure, which already develops secondary (counter) rolls in the corners, enjoys 
top-bottom symmetry. 



secondary roles are of considerable size for the BL separations of the up and down going 
flows, which approach the top and bottom plates, respectively, but are nearly invisible 
for the horizontal ones, which approach the side walls (at least for these Rayleigh num- 
bers). This might be attributed to the plume creations in the BLs on the bottom and 
top plates, which affect the horizontal but not the vertical sections of the flow. 

We emphasize that the convergence of the statistics implies that a center-point sym- 
metry should be established for the conditionally time averaged field u. Indeed, within 
5% precison we achieve such a center-point symmetry in our numerical simulations. As 
we shall see later this accuracy is by far sufficient to discriminate the main NOB effects 
in comparison to the OB results. 

For information we add the probability density function (PDF) of the center vorticity 
lu c , (nondimensionalized by the molecular vorticity v m L~ 2 ), see Fig 112] Already for Ra = 
10 6 two preferred values of the LSC can be recognized, reflecting clockwise and counter- 
clockwise rotation of the large scale convection roll. For Ra = 10 8 the two preferred 
vorticities are even more pronounced, reflected by the sharp peaks in the PDF. While in 
the chaotic phase (Ra ~ 10 6 ) the PDF is still broad, the flow cessations at Ra « 10 s are 
very fast events, since the small probability for ui c = suggests that the flow changes its 
direction more or less momentarily, leaving only a very small probability to find lo c = 0. 

5.3. Wind profiles: amplitudes 

Because even the conditionally time-averaged velocity field as shown in Fig. [T^] has a rich 
spatial structure, it is not immediately apparent how to define "a characteristic single 
wind velocity amplitude U" unambiguously. Several definitions have been proposed in 
the literature. In the present section we introduce various possible definitions for wind 
amplitudes U and corresponding Reynolds numbers Re = UL/v m , by investigating global 
measures as well as local ones, and also certain profiles. 

i. Global. A global measure for the strength (amplitude) of the convection can be based 
on the volume average. While the velocity average over the full volume is of course zero, 
we can consider the velocity rms average and obtain the energy based wind amplitude 

i" = ^« + «s) yj = {{ I M ) v M 

Apparently, both the primary center roll as well as the secondary rolls contribute to the 
value of U E . While the primary roll covers the whole interior, which is expected to be 
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Figure 13. Probability density function of the vorticity u> c at the cell center at Ra — 10 6 , 10 s . 

mixed by turbulence and therefore has essentially uniform temperature T c , the secondary 
rolls experience cither the cooler top or the hotter bottom regions only. Since under NOB 
conditions the viscosity deviations are different near the top and bottom (and, of course, 
from the bulk), the secondary rolls are expected to have different properties among each 
other as well as relative to the primary center roll. U E represents a well defined mixture 
of all of them. 

ii. Local. To deal with the top and bottom differences, local wind amplitudes may be 
introduced by considering the values of the conditionally time averaged field u(x) at 
particular spatial positions Xj in the flow field. Particular positions are e.g. those, where 
the conditionally time-averaged velocity field u has peak values TJ Pj along some vertical 
lines, labelled by j. In the following we shall consider such local peak amplitudes U Pj 
at peak positions on vertical lines with abscissas Xj = , -§ , and -| . These points Pj 
roughly correspond to the positions where the main, primary circulation roll is strongest 
(xi = ^p), to the flow maxima along the center line of the container (#2 = ■§), and to 
the region where the counter-rotating roll is well developed (x^ = j). 

iii. Profiles. Mixed type wind amplitudes, neither fully global nor fully local, can be 
introduced as area averages or rather, in 2D, as line averages. The areas (lines) can either 
be chosen as top/bottom plate parallel or as side wall parallel. The corresponding area 
averaged wind amplitudes then either, as in the case of plate parallel averaging, depend 
on the height z of the area (line) and lead to z-dependent wind profiles; this is relevant 
for the horizontal wind and its vertical profile U x {z). Or they depend, as in the case 
of vertical line averaging, on the x-distance to the side walls and lead to x-dependent 
profiles; these are relevant for the up rising or down falling flow and its horizontal profile, 
denoted as U z (x). 

Such area (line) and time averages are the relevant quantities in the well known re- 
lations between the dissipation rates e u , eg, and the nondimensionalized heat current 
density, i.e., the Nusselt number Nu. While the dissipation rates arc volume averages, 
the Nusselt number is defined in terms of a horizontal area average, which by conserva- 
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Figure 14. The vertical and horizontal velocity profiles U x (z) and U z {x) as derived from 
the conditionally time averaged velocity components and at three different Ra numbers 
10 4 ,10 6 ,10 8 . The abscissa and ordinate scales are the dimensionless width and height. The 
upper scales show the horizontal velocity profiles nondimensionalized with v m /L, i.e., Re x (z); 
the right scales show the vertical velocity's U z (x) profiles, also nondimensionalized by P m /L. 
Note the increase of the Re x<z scales with increasing Ra. Dashed lines indicate the OB case, full 
lines NOB case. In both cases A = 40K. Also indicated are the corresponding thermal slope BL 
widths Ajj t , which strongly decrease with Ra. 



tion of energy even is independent of height z. One contribution to Nu comes from the 
(Tu z ) J 4 j i-corrclation, the other one from the gradient of the temperature profile (T)A,t(z)- 
This justifies to introduce the mentioned area (line) averages also for velocity compo- 
nents. In particular, the vertical profile of the horizontal velocity is of interest, U x (z), as 
well as the horizontal profile of the vertical velocity, U z {x), precisely defined as 

U x {z) = {u^) x(z) 0), and U z (x) = {uZ) z{x) (x). (5.4) 

where (. . -)x(z) an d (■ ■ ■) z(x) represent the line averaging along the x direction for fixed 
z or along the z direction for fixed x. 

Area averaged profiles are displayed in figurc [T4l Surprisingly there are ranges with neg- 
ative (positive) U x (z) of the area averaged horizontal velocity in the immediate vicinity 
of the bottom (top) plates. Formally their origin is that the area averaged profiles U x (z) 
take notice of the sign of the corresponding local velocities and these in the secondary 
rolls are opposite to the center roll. Thus one clearly sees the effects of the secondary 
rolls in the inversion of the vertical profile in the neighborhood of the bottom and top 
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Figure 15. The vertical and horizontal rms velocity profiles t/£ ms (z) and Ul ma {x) as derived 
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plates. Physically this means that plumes in this range of small distances from the plates 
are mostly advected in the opposite direction until they come farther away. One also 
observes broken top-down symmetry. This is caused by the secondary rolls, since these 
cover regions of different temperatures. 

iv. rms profiles. In addition to area averages of the (conditionally time averaged) ve- 
locity components themselves one might also wish to analyze area averages of the rms 
fields of the corresponding components, defined as 

UT s {z) = ^(^) x{z) (*) = ^>, Wlt (*). and Ul™{x) = Jffi^ (*)• (5.5) 

These arc displayed in figure fT5l 

In the rms profiles U™ s one also observes the kink and the broken symmetry as 
in the averaged velocity profiles U x . z in figure [TH but no change in sign, of course. 
The preferentially strong plume convection near the side walls is clearly reflected in 
the U^ ms (x) profiles and their sharper peaks as well as apparently significantly smaller 
boundary layers than in the U^ ms (z) profiles. In the vertical profiles of U^ ms (z) the near- 
wall inversion (owing to the opposite sign of the velocity in the secondary and center rolls) 




is not present, in contrast to the amplitude averaged U x (z) profile in figurc[14l since U x ms 
is positive everywhere. Instead, one recognizes the signature of the secondary roll effect 
as the steep slope in the U" ns profile between the plate and its maximum position. 

Experimentally often Lagrangian flow properties as e.g. the plume turnover times are 
used to characterize in particular the large scale coherent flow (LSC) . We emphasize that 
such features should be identified in the U x (z), U z (x) profiles. The U x ms (z), U" as {x) pro- 
files, instead, reflect the energy strengths of the considered components. In the following 
wc will use two different amplitude velocities defined on the local maxima of horizontal 
profiles U x (z) and U x ms (z), denoted respectively as U^ 2 , U x ms,Mrms , where the super- 
script Mi indicates the the second (and positive) maximum on U x (z). In Figures [Til and 
[15] the z-distances, at which the velocities have the values U^ 2 and U™ s ^ M -r ms j respec- 
tively, are comparable. Both are located within the thermal BL at Ra = 10 4 while they 
are outside at Ra = 10 6 and 10 s . Furthermore, NOB effects on the maximum positions 
arc Rayleigh number independent, in the sense that one observes shifts of comparable size 
towards the bottom plate at Ra = 10 4 and 10 6 . However, as we shall see later, one no- 
tices slight differences because the nontrivial spatio-temporal flow structure is differently 
reflected in the U x and U x ms profiles. 

5.4. Scaling of amplitudes with Ra, Oberbeck-Boussinesq case 

Having described the flow structures and the definitions of several relevant measures for 
the magnitude of the thermally driven convection, we now offer our results on the Ra- 
and A-dcpendcncc of the various U amplitudes and the corresponding Reynolds numbers 
Re = U/{y m L~ 1 '). We start with the OB case, i.e., having temperature independent 
material parameters throughout the container, their values taken at the given arithmetic 
mean temperature T m . FigfTBI shows the scaling of Re E , Re 2 and Re rms - Mrms with Ra 
for the OB case. 

The structure in the Req B versus Ra curve around Ra « 10 5 is due to changes in 
the still present coherent flow patterns, implied by the boundary conditions. They are 
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10 4 10 s 10 6 10 7 10 ! 

Ra 



Figure 17. Reynolds numbers Re Pj at xj — (triangles), ^(crosses), (circles) in the OB 
case versus Ra for water at T m = 40° C; again A = 40K. The inset shows a sketch of the flow 
with the positions of the velocity peaks being indicated by the three different symbols. 



observed in the Nu behavi o r too (not shown here), and are also detected in experiment, 

" (|l975h . 



see for instance iThrelfall I (|1975f ). For water as the working fluid the typical spatial 
coherence length (in terms of L) of the coherent flow struct ures has decreased to a bout 
order 0.1 only in the i?a-range between some 10 7 and 10 8 ( Sugivama et 
power-law fit, Re E = cRa 1 , in the range 7 • 10 5 < Ra < 10 8 , corresponding to the 
transition range from chaotic to turbulent behavior, gives the exponent 7 = 0.616. 

The behavior of the Re Pj at Xj = j , |f , ^¥ ) see Fig- H3 is more noisy but scaling- 
wise similar to (at least compatible with) Re , although in Re Pj the spatial structures 
of the flow field are well taken into account. Re E is more robust with respect to the 
convergence of the statistics. Experim entally at moderate Ra the Reynolds numbers may 
scale differently, cf . lLam et al. ( 2002 ) . Here the Reynolds number based on the maximum 
horizontal velocity near the bottom plate outside the BL is reported to scale cx Ra 70 
for Ra < 2 • 10 7 and oc Ra 0A95 for Ra > 2 • 10 7 . The set of exponents 7 obtained in our 
simulation for Re E and the three Re P . (7 ~ 0.62) are therefore reasonably consistent 
with experimental findings. 

5.5. Scaling of amplitudes with Ra and A, non-Oberbeck-Boussinesq case 

We now study the NOB case, starting with the conditionally time averaged velocity field 

u, which now is bottom-top asymmetric. In particular, there is thinning/thickening of 

the bottom/top kinetic BLs. For the wind amplitudes based on the peak values of u along 

p p 

vertical lines we can distinguish between bottom and top peak velocities U b 3 and U t 3 



for the various x,-lines. As shown in Fig.[THl the ratio Re 



Pr=, 



I Re 



P~— T 



characteristic 



for the peak velocities taken in the main primary circulation roll converges to unity as 
Ra is increased up to Ra ss 10 8 and A kept fixed. This shows that for the primary roll, 
which scans the bulk of the convection cell, a single velocity amplitude develops in the 
tu rbulent range also under N OB conditions, just as assumed in the BL theory introduced 
in lAhlers et all ( 2006L 12007 ). This is also found in Figure HU where the peak velocity 



scales are replace by U x 1,2 and jj™s,M rm3 ^ j n particular, the large variation in the 
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Figure 18. Ratio -Re (r J JVOS /i?e t /J vos vs. Ra for the different peak positions in the cell. Note that 
due to the center-point symmetry of the conditionally time averaged field u, the peak positions 
along the lines Xj,b = ■§■,§> ^ for the bottom part have to be compared with the peak positions 
along Xj tt = \ for the top part. The figure shows that the main convection loop for large 

Ra establishes a uniform velocity amplitude, while the secondary counter-rotating rolls do not 
enjoy the same property. Instead, since they are BL-dominated rather than bulk-dominated, 
they show significant NOB bottom-top asymmetry. Note that the accuracy of this test may be 
assessed from the level of center-point asymmetry we had in the OB case, which was always 
below 5 %. Therefore we may conclude that in the chaotic regime there are no NOB deviations 
distinguishable for the Re b . 1 NOB /Re t . : ' NOB ratios at the two positions £2,3 = But there is 

a significant deviation for L/8, i.e., for the secondary rolls' ratio. 

bottom/top Reynolds number ratio in Figurc[THl which reflects the spatial inhomogeneity 
in the velocity distribution, is smoothened by taking the area average as indicated by 
the almost flat profiles of U* u and U r x ms ' M ^° in Figure [TO] 

The secondary rolls, counter-rotating to the primary roll, have different wind ampli- 
tudes near bottom and top, as becomes apparent in the bottom-top ratio at the line 
Xib/L — 0.125 and its mirror at the top Xi^/L = 1 — 0.125. The secondary roll at the 
warmer bottom plate is faster than the corresponding secondary roll near the top plate. 
We understand this from the smaller viscosity near the bottom due to the higher temper- 
ature. For the working fluid water and A = 40K the bottom-top asymmetry is as large 
as w 25% 

We therefore calculated the NOB/OB ratio of Re E versus A for various values of Ra 
as well as versus Ra at fixed non-Boussineqness A = 40K. This is shown in Fig. [201 We 
observe that NOB effects arc clearly present, although rather weak only, of order 2%, 
indicating a small increase of the kinetic energy based mean velocity. The reason for the 
rather small NOB effect on the global wind amplitude in spite of the large changes of the 
bottom and top velocities is that the secondary rolls only contribute a limited fraction 
to the global volume average. Also, the changes of the velocities of the secondary rolls 
are opposite in sign, one contributing a larger (bottom) the other contributing a smaller 
(top) amplitude. The remaining net change of the global amplitude U E thus again is 
due to the nonlinear temperature dependences of the material parameters, producing 
different secondary roll velocities at bottom and top. This crucial importance of the 
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Figure 19. Same as Figure [TH] but with the peak values of the horizontal area averaged 



velocities lit' 2 and UZ 



taken for the velocity scales. 



nonlinca rities in the t empe rature dependence of the material properties was already 
found in lAhlers eiaJ] (|2006l) . 



Let us recall that the Nusselt number becomes smaller with increasing A as was shown 
in Fig. [SI Therefore, although the global Reynolds number ratio is enhanced under NOB 
conditions, the overall heat transport is attenuated. We have also included in Fig. [20] a 
comparison with BL theory. We remind that to make this theory predictive an additional 
assumption on F\ (e.g. F\ ~ 1) has to be made. Note that the extended BL theory in 
this form undcr-estimates the NOB effect on the wind. 

As the origin of this discrepancy we can now identify the effect of the temperature 
depend ence of the thermal expansion coefficient, which is not included in the BL equa- 
tions in Ahlers et al. ( 20061 ). For this explanation we offer the following argument. As a 



naive estimate one can assume that the volume averaged velocity scale U E in essence 
should coincide with the free-fall velocity, that is U E ~ ^ gLf3(T)A. This gives the scal- 
ing Re ~ Pr^ 1 / 2 Ra 1 / 2 . Assuming then that the temperature of the bulk is dominated 
by T m and T c respectively in the OB and NOB cases, one gets 



E /fl/TU V2 

(5.6) 



RtNOB ( /TO 



This is in encouraging agreement with the DNS data, cf. Figurc[20l - As a further support 
of our argument emphasizing the importance of the temperature dependence of the ther- 
mal expansion coefficient in the bulk we also calculated Reff OB / Re E B for a hypothetical 
liquid, which has all material properties as water, apart from the thermal expansion co- 
efficient (3, which we keep constant at /3 m , see Fig. [21] Indeed, Re E TOB now only shows 
a smaller than 1% deviation from ReQ B , even at A = 60K. These tiny deviations from 
R c nob/ R e OB = 1 are consistent with the results from the Prandtl Blasius theory with 
the additional assumption F\ = 1. Fig.[2"T1thus confirms that the main origin of the NOB 
deviation in the Reynolds number Re E is the temperature dependence of the thermal 
expansion coefficient, an effect which clearly cannot be captured in the Prandtl-Blasius 
BL theory. 
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Figure 20 
water at T, 



The Reynolds number ratio RefiQB / R e OB based on the total kinetic energy for 
= 40°C. Upper: Ref 



Re NOB / Re OB versus Ra at fixed non-Oberbeck-Boussinesqness A 



/Re B versus A for various choices of fixed Ra. Lower: 

40K. The full line is calcu- 
lated from BL theory under the assumption that Xf, +\f = 2\q B , i.e., Fa = 1. The dashed line 
indicates the limit A — > where the Reynolds number ratio approaches 1. For the explanation 
in terms of the T-dependence of j3 see text, in particular (|5.6|l . 



The changes of Re caused by the loss of OB conditions can be analysed quantitatively 
in detail. For this we decompose the volume average into three regions, corresponding to 
the main, primary, large scale circulation (LSC) and the two secondary lower and upper 
corner rolls. 

(Re E ) 2 = a LS c(Ref sc ) 2 + ai ower (Ref ower ) 2 + a upper (Re^ pper ) 2 . (5.7) 

Here a a denotes the ratio of the volume V a occupied by the region a, with a = LSC, 
lower, or upper secondary roll, to the total volume, and Re E is the Reynolds number 
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Figure 21. Same as Figure l20l upper, but here with only linear dependence of p(T) on the 
temperature T, meaning that the thermal expansion coefficient (3(T) is set to the constant value 
(3 m . Note that under this approximation there is good agreement with the results from BL 
theory and F\ = 1 for large Ra, which means that the thermal convection takes significant 
notice of the full T-dependence of the thermal expansion coefficient f). BL theory misses that 
per construction. 



based on the kinetic energy averaged over each volume V a , 




Re£ = ^l with U^J-l^x-{ul(x)+ul(x)). 

The lower secondary roll is characterized by a clockwise rotation and negative values of 
the tZ(a;)-stream function ^{x, z) = J Q dz u^{x, z). The region Vi ower of the lower corner 
roll is defined by tp < 0, x < L/2 and z < L/2, while for the upper corner roll we have 
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Figure 22. The relative deviation of the kinetic energy {a a ,NOB{R&a,NOBY !_a a,os(-Ref ob) 2 } 
/(Reo B ) 2 in region a due to NOB effects for water at Ra = 10 8 and T m = 40°C. The label a 
denotes the portion of the large scale circulation (LSC) roll in the center region (as indicated 
by the solid lines in Figure [12]) . or the secondary counter-rotating rolls in the lower and upper 
corners (given by the dashed lines in Figure [12}. The partition is determined by using the sign of 
%j) = f* dz u^ c {x, z), e.g., the lower corner flow is defined as the region satisfying ip < 0, x < L/2 
and z < L/2. Upper panel: the full temperature dependence is considered for the buoyancy 
g(l — p/pm) as given in Table [1] Lower panel: restriction to only linear dependence of the 
buoyancy with respect to the temperature T, i.e., the thermal expansion coefficient /3(T) = f3 m 
is constant. The insets show the (squared) NOB/OB Reynolds number ratio in the each region. 



ip < 0, x > L/2 and z > L/2. The remaining region comprises the primary main roll, 
LSC. 

The insets of Figurel2"2"lshow the (squared) NOB/OB Reynolds number ratios (Re-^ nob) 2 
I (R e a ob) 2 f° r eacn region a at Ra = 10 8 . The insets both in the upper and lower pan- 
els reveal that the largest NOB enhancement of the kinetic energy occurs inside the 
lower corner secondary roll (near the warmer bottom plate), as already found in the 
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bottom-top asymmetry of the peak velocity Re^NOi! /^ e tTOB 8 m Figure [18] But even 
this enhancement does not impact significantly on the overall i?e s -change because the 
volume ratio of each corner flow is only a\ ower = a upper « 9% while the main vol- 
ume fraction is cllsc ~ 82%. This is even more pronounced in the OB case, where 
ai OW er{Ref ower ) 2 /{Re E ) 2 = 0.05 and a LS c(Ref sc ) 2 / (Re E ) 2 = 0.90, while for the NOB 
case at A = 40K these fractions are 0.07 and 0.04 for the lower and upper secondary rolls, 
and the primary LSC contributes 0.89. To visualize these contributions of each subvolumc 
a to the overall change, we show in Figure [22] the normalized kinetic energy deviation 
{aNOB, a {Refj OBjQ ) 2 - a B,a{ReoB. a ) 2 }l{ Re OB) 2 due to tlie N0B effect. The upper 
panel proves that the NOB enhancement of the total kinetic energy is primarily due to 
the LSC and secondarily to the lower corner roll. For comparison we plotted the corre- 
sponding deviations in the lower panel for a hypothetical fluid with thermal expansion 
coefficient fixed at (3 m . Then the enhancement of the LSC contribution is much smaller, 
while the attenuation and enhancement respectively in the upper and lower secondary 
rolls are comparable and in addition compensate each other. This leads to the much 
smaller change in the total kinetic energy in Figure [5T] as compared with allowing the 
full temperature dependence of /3(T) in Figure [2"D"1 Apparently the nonlinear temperature 
dependence of the buoyancy is very important. Note that the derivative of the buoyancy 
d(g(l — p/pm)) /dT, corresponding to the driving force per temperature displacement, 
increases when increasing the temperature deviation (T— T m ) because the coefficient C2 
for the buoyancy expression as reported in Table [T] is positive. Therefore the buoyancy 
force gets larger and the bulk kinetic energy is more enhanced for given temperature de- 
viation T — T c , if the mean bulk temperature T c is larger than T m as observed in Figures 
OandU 



6. Summary and conclusions 

In summary, we have studied the temperature profile, the heat current density, and 
the properties of the large scale convections as defined by several representative velocity 
scales. The center temperature T c and the Nusselt number ratio Nunob/Nuob result- 
ing from the two-dimensional numerical N OB simulations are in good agreem e nt wi th 
the available experimental data for water ( Ahlers et al. ( 20061) ). lAhlers et al. I (l2006lV s 
experimental finding F\ = 1 for water is argued to be incidental, originating from the 
specific temperature dependence of the material constants of water at 40°C. This finding 
cannot be generalized to other fluids or to other mean temperatures. For water the heat 
flux reduction due to the deviations from OB conditions is for all practical purposes due 
to the modified temperature drops over the BLs, represented by F&, whereas for other 
working fluids it is influenced also by the changes of thermal BL thicknesses, expressed 
by F x . 



The results of the simulations also agree with Ahlers et al. ( 20061 ) 's predictions for 



the central temperature, which is based on an extended Prandtl-Blasius theory. As that 
theory ignores plumes and side wall effects, these apparently hardly contribute to the 
determination of the central temperature T c . However, physically they do contribute to 
the shape of the temperature profiles. Our simulations reveal their increasing effect on 
the profiles with increasing Ra. As the overall heat transfer is determined by the slope 
of the temperature profiles at the plates, it is to be expected that w ith increasing Ra the 
plumes increasingly affe ct the Nusselt number, in coherence with iGrossmann fc Lohse 



(|2000l l200ll . 120021 (20041 ) 's unifying theory. 



The second part of the paper is devoted to the flow organization in the OB and the NOB 
cases. First of all, also in the NOB case the large scale convection roll is characterized 
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by only one velocity scale. In contrast, the top and bottom corner flows have different 
velocity scales in the NOB case, reflecting the enhanced and reduced viscosities close to 
the respective plates. We defined various different velocity scales, based on global and area 
averages and peaks in the profiles and analyzed how these change under NOB conditions. 
For the ratio of the energy based Reynolds numbers which also is representative for the 
others we find Re% OB /Re% B « (/3(T c )//3(T m )) 1 / 2 , i.e., NOB deviations in the Reynolds 
number are strongly governed by the temperature dependence of the thermal expansion 
coefficient. This finding suggests that fluids which display no or only a weak temperature 
dependence of the thermal expansion coefficient hardly show any NOB effects on the 
Reynolds numbers. 

From our point of view the two next steps in numerical work on NOB correction 
for RB flow arc: (i) Study the detailed modifications of the various BL thicknesses and 
profiles through NOB effects, and (ii) confirm that at least for Prandtl number of 1 
and larger the findings of this paper also hold for three-dimensional RB flow. Moreover, 
NOB experiments focusing on the flow organization, BL layers, and Reynolds number 
modifications would be very desirable. 
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